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Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are complied from the 

literature. Multiple optical emission components are extracted with power-law function 

CN ' fits to these lightcurves. We present a systematical analysis for statistical properties and 

^J their relations to prompt gamma-ray emission and X-ray afterglow for each component. 

We show that peak luminosity in the prompt and late flares are correlated and the 
evolution of the peak luminosity may signal the evolution of the accretion rate. No tight 
correlation between the shallow decay phase /plateau and prompt gamma-ray emission is 
found. Assuming that they are due to a long-lasting wind injected by a compact object, 
we show that the injected behavior favors the scenarios of a long-lasting wind powered by 
a Poynting flux from a black hole via the Blandford-Znajek mechanism fed by fall-back 
mass or by the spin-down energy release of a magnetar after the main burst episode. The 



as peaking later. Assuming that the onset bump is due to the fireball deceleration by 
the external medium, we examine the To — E~ » S o relation and find that it is confirmed 
with the current sample. Optical re-brightening is observed in 30 GRBs in our sample. 
It shares the same relation between the width and the peak time as found in the onset 
bump, but no clear correlation between L_r iP and E~^ so similar to that observed for the 
onset bumps is found. Although its peak luminosity also decays with time, the slope is 
much shallower than that of the onset peak, as is the case for the onset bumps. We get 
L oc £p , being consistent with off-axis observations to an expanding external fireball 
in a wind-like circum medium. Therefore, the late re-brightening may signal another jet 
component. Mixing of different emission components may be the reason for the observed 
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chromatic breaks of the shallow decay segment in different energy bands. 
Keywords: gamma-rays: bursts - methods: statistics 

1. Introduction 

Gamma-ray bursts (GRBs) and their broadband afterglows are the most luminous 
phenomenon in the Universe. The most popular model is the internal + external 
shock fireball model, which suggests that the prompt emission is produced by in- 
ternal shocks at a distance internal to the deceleration radius of the GRB fireball 
and the broadband afterglows are from the external shocks when the fireball is 
decelerated by the ambient mediu m I ^ Prompt and afterglow emissions involve 
two distinct processes at different sites in this model. The observations with Swift 
mission significantly improved our understanding of the internal+external shock 
picture for GRBs^H Rapid localization with Swift also revolutionized the ground- 
based follow-up observations, leading to establishing a large sample of GRBs with 
well-sampled optical afterglow lightcurve s and red shift measurement. Following our 
comprehensive analysis of the Swift dat a ^" I ' l " 1 ^ 1, we make a systematical analysis 
of the optical data and explore their relations to the prompt gamma-ray emission 
and the X-ray afterglow. 

2. Data 

We include all the GRBs that have optical afterglow detection in our sample. A 
sample of 225 optical lightcurves are complied from the literature in the period from 
Feb. 28, 1997 to November 2011. We make an extensive search for the optical data 
from published papers or from GCN Circulars in case of no published paper being 
available for some GRBs. Well-sampled lightcurves are available for 146 GRBs. We 
collect the optical spectral index (/3ojfJ and the extinction Ay by the host galaxy 
of each burst from the same literature. Galactic extinction correction is made by 
using a reddening map presented by Schlegel et al. (1998) . Since the Ay values 
are available only for some GRBs and the Ay is derived from the spectral fits by 
using different extinction curves, we do not make correction for the extinction of the 
GRB host galaxy. The fc-correction in magnitude is calculated by k = — 2.5(/?o ~ 
l)log(l + z). Note that most of the well-sampled optical lightcurves are in the R 
band. For a few GRBs, the data are well-sampled in the other bands. We correct 
these lightcurves to the R band with the optical spectral indices. For data in late 
epochs (~ 10 6 seconds after the GRB trigger), possible flux contribution from the 
host galaxy is also subtracted. The isotropic gamma-ray energy (-E^iso) is derived 
in the energy band of 1 — 10 4 keV in the burst local frame with spectral indices. 
Data and references of our full sample will be reported in a series of papers (in 
preparation). 



a An optical spectral index /3q = 0.75 is used for those GRBs without /3q available 
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3. Lightcurve Fitting and a Synthesized Optical Emission 
Lightcurve 

The optical lightcurves are usually composed of one or several power-law segments 
as well as flares/re-brightening features. The mix of different components makes 
the diversity of the optical afterglow light curv es. Different from previous statistical 
analysis on the optical data by some team a ' ' ^ ' ' ^ we fit the lightcurves with 
a model of several components in order to subtract each emission component from 
the lightcurves. The basic components in our model are a power-law function and 
a smooth broken power-law, i.e., 

F = F t' a , F = F Q [(t/t b ) a ^ + (t/t b ) a n~ 1/ul , (1) 

where a is the temporal decay slope, t b is the break (or the peak) time, and uj 
measures the sharpness of a break (or a peak) of a lightcurve. The width of a 
flare/bump is measured with the full-width-at-half-maximum (FWHM). We develop 
an IDL code to make the best fit with a subroutine called MPFITO. Note that the 
parameter lu is usually fixed at 3 or 1 in our fitting. The approach of our lightcurve 
fitting is as follows. Initially, we add components to our model by inspecting the 
global feature of a lightcurves. If the reduced x? is much larger than 1, we continue 
to add components and make the fit. We try to get a fit with \\ being close to 1. The 
Xr values for some lightcurves are much lower than 1, indicating that some model 
parameters are poorly constrained. Therefore, we fix some parameters to make the 
fits for these GRBs. The erratic fluctuation of some data points with a small error 
bar in some GRBs, such as GRB 030329, makes x^ De much larger than 1. We do 
not add additional component for these data points and the x? of the fits for these 
GRBs are much larger than 1. The most difficult problem of our fit is extraction 
of the seriously overlapped flares/bumps from the lightcurves. The slopes of these 
flares/bumps are usually quite uncertain. In our fitting, we first let all parameters 
free to get the best fit for the global lightcurve, then adjust the rising slopes to 
ensure that the fitting curve crosses the data point around the peak times of the 
two components. Finally, we fix the rising slopes and make the best fit again. 

The flux of an flare event usually rapidly increases and drops. We identify a flare 
event with a criterion that the slopes of the rising and decaying parts are steeper 
than 2. Those optical flares during the burst duration are defined as prompt optical 
flares. Reversed shock flares following the prompt flares are observed in only a few 
GRBs. They are extensively discussed in literature. The standard fireball model 
suggests that the decay slope of the afterglows should be steeper than 0.75 if no any 
late energy injection after the GRB phase. We therefore define the shallow decay 
phase with the criterion that the initial decay slope of this segment is shallower than 
0.75, which transits to a steeper decay after the break time. An optical afterglow 
onset is defined as an initial smooth hump with a peak less than 1 hour after 



b http://www. physics. wisc.edu/ craigm/idl/ntting.html 
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the GRB trigger and decaying as a power-law with a slope being consistent with 
the external shock models (0.75 < a < 2). A re-brightening hump is analogous 
to the afterglow onset hump but is later than the onset hump. The supernovae 
bumps are identified as those late re-brightenings that peak at around 1 ~ 2 weeks 
post the GRB trigger. We decompose emission components and make statistical 
analyses for each one. Morphologically, a synthesized optical lightcurve of these 
components along with some typical examples is shown in a cartoon picture (Figure 
[1]) based our statistical results. It describes 7 components related to the X-ray 
canonical lightcurve^. The early optical afterglow lightcurves (t < 10 3 second post 
the GRB trigger) of about one-third GRBs show a smooth bump and the other one- 
third lightcurves start with a shallow decay segment. Twenty- four optical flares are 
observed in 19 GRBs. Late re-brightenning is observed in 30 GRBs. A jet like break, 
in which the decay slope transits from 0.75 ~ 1.5 to 1.5 ~ 2.5, is detected in 10 
GRBsjj. A clear Supernova bump is detected for 18 GRBs. The detection probability 
of each component is also marked in the cartoon lightcurve. We report our results 
for the flares, shallow decays, afterglow onset bumps, and late re-brightening in 
this paper. We mark the parameters of these components with superscripts "F" 
for the flares, "S" for the shallow decays, "A" for the afterglow, and "R" for the 
re-brightening. 



b. Reverse Shock 
.-.(10/146) IV. Late Optical Flare 
.'', (19/146) 



'. II. Shallow Decay 
j (39/146) 

., V. Re-Brightenir 

(30/146) 




Time since trigger(s) 

(a) A synthesized cartoon lightcurve of mul- 
tiple optical emission components based on 
our statistics, (b) Examples of the lightcurves 
with various emission components. The solid 
lines represent the best fit to the data. 
Simultaneous X-ray data observed with 
Swift/X.HT (crosses with error bars) are also 
presented. 
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Fig. 1. 



(b) 



c We do not include those breaks that have a slope shallower than 0.75 breaking to a slope steepen 
than 1.5 
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4. Flares 

We get 24 late flares in 19 GRBs. Relations of the width (w F ) and the peak lumi- 
nosity (L F iso ) of the flares as a function of the peak time (i_) are shown in Figure 
[21 The t F ranges from ~ tens of seconds to ~ 10 6 seconds. The w F values are in 
the same range as t F . The L^ iso ranges in 10 43 — 10 49 erg s _1 , with a typical value 
of 10 46 erg s _1 . A tight correlation between w F and t F is found. The best fit gives 
logw F = -0.32 + 1.01 log i^, i.e., w F ~ t%/2. The L F Rp is anti-correlated with t F 
in the burst frame, i.e., logL^ iso48 = (1.89 ± 0.52) - (1.15 ± 0.15) log[t£/(l + z)] 
with a Spearman correlation coefficient of 0.85 and a chance probability p < 10~ 4 . 
Therefore, a flare peaking at a later time tends to be dimmer and wider. 

-^r iso are usually smaller than 1/100 of S 7i ; so . The L^ iso is correlated with 
L 7 , iso ,'i.e., logLJ^/10 48 = (-3.97 ±0.60) + (1.14 ± 0.27) log L 7 , iso /10 50 with a 
Spearman correlation coefficient of r = 0.75 and a chance probability p ~ 10~ 3 . The 
flares in GRBs 050401, 060926, and 090726 are out of the 3ct region of the fit. With- 
out considering the flares in the three GRBs, it is found that the t F is also tightly 
anti-correlated with E lMo i.e., logt'J = (5.38 ± 0.30) - (0.78 ± 0.09) \ogE lMo /10 5Q 
(with r = 0.92). Similarly, a tight anti-correlation between £r, p and t p in the 
burst frame is found without considering the flares in the three GRBs, e.g., 
log[ip7(l + z)] = (7.57 ± 0.60) - (1.35 ± 0.17) log £ 7 ,i SO ,50 with a Spearman cor- 
relation coefficient of 0.91. 

It is interesting to study whether the optical flares are associated with X-ray 
flares. Early optical flares are only observed in the lightcurves of GRBs 060210, 
060926, 090618, and 090726 in our sample, indicating that the fraction of GRBs 
with detection of early optical flares is much lower than that of the X-ray flares. 
Among the 19 GRBs with detection of the optical flares 16 are detected with Swift. 
Simultaneous observations with XRT during the optical flares are available for 
GRBs 050401, 060206, 060210, 060607A, 060926, 070311, 071010A, 071031, 080506, 
090618, and 100728B. An X-ray flare that may be associated with the optical flare is 
only observed in GRBs 060926, 070311, and 071010A. The optical flares of the three 
GRBs are lagged behind the corresponding X-ray flares. Measuring the lags with 
the peak times of the flares, we get 196 seconds, 7.7 x 10 4 seconds, and 2.45 x 10 4 
seconds for the flares in GRBs 060926, 070311, and 071010A, respectively. The lag 
is potentially proportional to the peak time of the flares with the three flares. 



4.1. Prompt Optical and Reversed Shock Flare 

Well-sampled prompt optical flares are observed for GRBs 061121, 060526, 080129, 
080319B, and 110215A. They usually trace the pulses of the prompt gamma-ray 
phase with a significant temporal lag. Reversed shock Flares are detected for GRBs 
990123 and 060111B. Generally speaking, the light curve dominated by reverse the 
shock emission is fast rising and decaying, quite similar to the prompt optical flares. 



February 26, 2013 20:51 WSPC/INSTRUCTION FILE MS 'Enwer Liang 



6 Liang et al. 
4.2. Late Flare 
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Fig. 2. Correlation between w F and K as well as relations of L^ iso to t£ and L 7i i so - The black 
solid dots, grey open circles, and grey triangles are for the optical flares, X-ray flares, and prompt 
gamma-ray pulses. Best fit lines with 3<r significance level are also shown. 



5. Early Shallow Decay Segment 

We get a sample of 41 GRBs with a shallow decay segment from the 146 GRBs. 
Thirty-one out of the 41 shallow decay segments transit to a decay slope of 1 ~ 2.5, 
and 5 of them are followed by a sharp drop with a decay slope being steeper than 
2.5. About half of the shallow decay segments look like a plateau, with \af J < 0.3. 
Figure [3] shows the distributions of the break times (i^) and the luminosity at 
the break (Lf iso ) as well as their correlation. The break time ranges from tens of 
seconds to several days post the GRB trigger, with a typical t^ of ~ 10 4 seconds. 
The L^ h varies from 10 43 to 10 47 erg s _1 , and even ~ 10 49 erg s _1 of the early 
break in some GRBs. The L^ b is anti-correlated with t^, which is logL| 48 = 
(1.75 ± 0.22) - (0.78 ± 0.08) log[^/(l + z)] with a Spearman correlation coefficient 
of r = 0.86 and p < 10" 4 . 

A shallow decay segment is commonly seen in the well-sampled of XRT 
lightcurves, except for a few GRBs that their XRT lightcurves decay as a single 
power-law-^ (Liang et al. 2009). It was also reported that the X-ray luminosity at 
the break time is correlated with the break time L^Dainotti et al. 2010). We over- 
plot the Lb, iso as a function of % in the burst frame in Figure El One can observe 
that optical data share the same relation to the X-ray data. Note that the X-ray 
luminosity is in the 0.3-10 KeV energy band and the optical data are in the R band, 
so that the X-ray luminosities are significantly higher than the optical ones. The 
observed photon indices of the X-ray spectra are ~ 2. Therefore, the energy spectra 
of the X-rays are flat and the derived L D ,iso — ^b relation in the 1 KeV band is 
roughly consistent with that observed in the X-ray band. 
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We examine the chromaticity of the shallow decay segments in the X-ray and 
optical bands. The X-ray observations are available for 17 out of the 34 GRBs. We 
extract the underlaying afterglow components from the X-ray data and compare 
the af , a§ , and t^ in the X-ray and optical bands in Figure [4] It is found that 
the data points are scattered around the equality line, and a tentative correlation 
between the break times of the optical and X-ray lightcurves is observed, with a 
chance probability of the correlation of ~ 0.15. These is no correlation between the 
decay slopes of the X-ray and optical lightcurves. The decay segment prior to the 
break times in the optical bands tend to be steeper than that of the X-ray band, 
but the slopes post the breaks are roughly consistent, except for those af > 2.5 in 
the optical bands. 
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Fig. 3. Correlation and Distributions of L^ iso and if /(l + z) of the shallow decay segments for 
the GRBs in our sample. The grey circles are for the the X-ray data from Dainotti et al. (2010). 
Lines arc the best fit line. 
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Fig. 4. Comparisons of the decay slopes and the break times in the optical and X-ray bands. The 
dashed lines are the equality lines. 
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6. Early Afterglow Onset Bump 

An early smooth bump is observed in the optical afterglow lightcurves of 42 GRBs 
in our sample. The peak luminosity and width as a function of t~ are shown in 
Figure [5] It is found that L^ is anti-correlated with tp measured in the burst 
frame and w is tight correlated with t^ , indicating that it is wider and dimmer if it 
peaks later. The isotropic prompt gamma-ray energy (i? 7j i so ) is tightly correlated 
with L^ p (Figure O. The best fit yields L^ lso oc £^?° o ±0 - 14 . Assuming that the 
bumps signal the deceleration of the GRB fireballs in a constant density medium, 
we calculate the initial Lorentz factor (Fo) of the GRBs with redshift measurements. 
The derived Tq are typically a few hundreds. The Tq — -E 7] ; so relation discovered by 
Liang et al. (2010/^ is confirmed with the current sample. 




10' 10' 10' 10* 10' 10 ! 10' 10' 
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Fig. 5. Correlation between w and t^ as well as relations of L^ iso to t^ and E 7i i ao . The black 
solid dots are for the afterglow onset bumps. The data for the late re-brightenings are also shown 
with opened triangles for comparison. 



7. Late Re-Brightening 

A re-brightening hump is analogous to the afterglow onset hump but it is latter 
than the onset humps. It is observed in 30 GRBs in our sample. Both a\ and OL2 of 
the onset and re-brightening bumps are well consistent. Most a^s are in the range of 
—3^0, with a typical value of ~ —1. The ot\ in some cases is smaller than —3. The 
oli values are in the range of 0.3 ~ 2. It is found that the £?■ randomly ranges from 
several hundreds of seconds to several days. The Lp is systematically lower than 
Lp . The typical E§ iso is 10 45 ~ 10 4r erg. Correlation of the characteristics of the 
re-brightening bumps are also shown in Figure [5] in comparison with the afterglow 
onset bumps. It shares the same relation between the width and the peak time as 
the onset bump, but no clear correlation between Lr^ p and E lt i SO , as is the case 
for the onset bumps. Although its peak luminosity also decays with time, the slope 
is much shallower than that of the onset peak. We get L ex t~ , being consistent 
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with off-axis observations of an expanding fireball in a wind-like circum medium 1 -^. 
Therefore, the late re-brightening may signal another jet component. 

8. Summary 

We have analyzed a well-sampled lightcurves of 146 GRBs. A generic optical 
lightcurve and statistical results of various emission components are presented. We 
summary our results below. 

(1) Twenty-four late optical flares are obtained from 19 GRBs. The fraction of 
the detected optical flares is much smaller than that of X-ray flares. Associated 
X-ray flares are observed for 4 optical flares and the optical flares usually lag be- 
hind the corresponding X-ray flares. We find L^ iso oc Ll:\l °' 27 , w F ~ ip/2 and 
^R.iso ^ [*p/(l + z)] _1 - 15±0 - 15 , indicating that the optical flares are correlated with 
the prompt gamma-ray phase and a flare peaking latter is wider and dimmer. These 
results suggest that the physical origin of the late optical flares could be the same 
as the prompt gamma-ray phase and the temporal evolution from the GRB phase 
to late optical flares may signal the global evolution of the GRB central engine. 

(2) A shallow decay segment is observed in 39 GRBs. The detection fraction 
of the optical shallow decay component is comparable to that in the X-ray band. 
The X-ray and optical breaks are usually chromatic, but a tentative correlation 
is found. Their break times (ijjj) range from tens of seconds to several days post 
the GRB trigger, with a typical value of ~ 10 4 seconds. The break luminosity is 
anti-correlated with t^, Lj| b oc [t^/(l + z)\ ' , similar to that derived from X-ray 
flares. The shallow decays / internal plateaus may be evidence of a long-lasting wind 
powered by the central engine. The injection behavior may be used to diagnose the 
nature of the central objects in the GRB central engine. Assuming that the behavior 
of the luminosity injected into the forward shocks evolves as L — Lot~ q , we find 
that the long-lasting wind may be powered by a Poynting flux from a black hole 
via the Blandford-Znajek mechanism fed by fall-back m ass or by th e spin-down 
energy release of a magnetar after the main burst episod e ll1 I One critical 
issue to explain the shallow decay segment with the energy injection scenario is 
the chromatic breaks in the optical and X-ray bands. Mixing of different emission 
components may be the reason for the observed chromatic breaks of the shallow 
decay segment in different energy bands. 

(3) An early smooth bump is observed in the optical afterglow lightcurves of 
42 GRBs in our sample. It is found that £r p is anti-correlated with t£ measured 
in the burst frame and w is tightly correlated with t^, indicating that a dimer 
flare tends to peak later and be wider. The 2? 7) is is also tightly correlated with 
L^ . Assuming that the bumps signal the deceleration of the GRB fireballs in a 
constant density medium, we calculate the initial Lorentz factor (Tq of the GRBs 
with redshift measurements. The derived Tq are typically a few hundreds. The 
To — -E 7 ,iso correlation discovered by Liang et al. (2010p^ is confirmed with the 
current sample. The tight relation of the onset bumps to the prompt gamma-rays 
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may open a window to investigate the radiation physics of GRB fireballs. 

(4)A re-brightening hump is analogous to the afterglow onset hump but it follows 
an onset hump or a power-law decay segment. It is observed in 30 GRBs in our 
sample. It shares the same relation between the width and the peak time as for the 
onset bumps, but no clear correlation between L^ p and E^^ so is found. Although 
its peak luminosity also decays with time, the slope is much shallower than that 
of the onset peak. We get L oc t" 1 , being consistent with off-axis observations to 
an expanding external fireball in a wind-like circum medium. Therefore, the late 
re- brightening may signal another jet component. 
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